%最小二乘拟合
%x,y为待拟合得数据，返回值p为拟合函数的系数（按x的升幂排列）
function [p]=LSM(x,y)
if length(x) ~= length(y)
warning('输入矩阵维度有误');
end
g0=x.^0; 
g1=x;
g2=x.^2;
G=[g0 g1 g2]; %G矩阵构造
%=====normal equation
A=G'*G; %法方程组的系数矩阵
b=G'*y;
p=A\b;  %求解法方程组
end
